Weak values are universal in von Neumann measurements 



Justin Dressel and Andrew N. Jordan 
Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, USA 

(Dated: March 8, 2013) 

We refute the widely held belief that the quantum weak value necessarily pertains to weak mea- 
surements. To accomplish this, we use the transverse position of a beam as the detector for the 
conditioned von Neumann measurement of a system observable. For any coupling strength, any 
initial states, and any choice of conditioning, the averages of the detector position and momentum 
are completely described by the real parts of three generalized weak values in the joint Hilbert space. 
Higher-order detector moments also have similar weak value expansions. Using the Wigner distri- 
bution of the initial detector state, we find compact expressions for these weak values within the 
reduced system Hilbert space. As an application of the approach, we show that for any Hermite- 
Gauss mode of a paraxial beam-like detector these expressions reduce to the real and imaginary 
parts of a single system weak value plus an additional weak-value-like contribution that only affects 
the momentum shift. 
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Since its introduction in 1988 by Aharonov, Albert, 
and Vaidman (AAV) [T] and subsequent confirmation 
[21 [3] , the weak value of a quantum observable has been 
a source of considerable controversy. AAV showed that 
a weak conditioned von Neumann measurement which 
coupled an observable A to a continuous detector con- 
sistently produced the complex weak value expression, 
(A) w = (i{;f\A\ipi) /(ipflipi) in the detector's linear re- 
sponse after pre-selecting the system state to \tpi) and 
post-selecting the system state to \tpf). Notably, the 
parts of this complex expression need not be constrained 
to the eigenvalue range of A, a fact which has prompted 
considerable recent interest both for amplifying the mea- 
surements of small quantities in weak measurements [ll [5] 
and for fruitfully using weak measurements to interpret 
quantum phenomena [5lH2| . 

There has also been considerable recent interest in gen- 
eralizing the derivation of pre- and post-selected mea- 
surements beyond the weak measurement regime consid- 
ered by AAV. Example efforts include the increase of the 
coupling strength [TUl EJ EH [H] , the addition of detector 
dynamics [TU [TS] , the addition of decoherence and noise 
[P7] . treatments of orthogonal post-selections [Hj, con- 
siderations of full counting statistics [l9j, a realization 
with Fock states [20], and the determination of optimal 
detector states [2"Tj . The AAV regime weak value has 
also been generalized to mixed initial states pi and ar- 
bitrary post-selections represented by positive operators 
Pf [22H24], 
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Notably, Eq. ([IJ reduces to the original expression when 
Pi = \ipi){ipi\ an d Pf — IV'/X^/L but also has the benefit 
of subsuming the expectation value of A as a special case 
when Pf = 1. 

In this Letter, we extend these works with five main 



results. Our primary result is to show that all von Neu- 
mann measurements are exactly described by generalized 
weak values such as Eq. for any coupling strength, any 
choice of initial mixed system or detector states, and any 
choice of generalized post-selection. Hence, weak values 
are universal in von Neumann measurements, and thus 
are not solely peculiarities of the AAV weak measure- 
ment regime. Our second and third results are compact 
expressions for the relevant generalized weak values in 
terms of the Wigner distribution of the detector. Fi- 
nally, our fourth and fifth results are applications of our 
general results to transverse Hermite-Gaussian modes of 
a detecting beam, such as those naturally produced by 
laser cavities. In the Supplementary Material [25] we fur- 
ther generalize our main results to higher-order detector 
moments and arbitrary Hermite-Gauss detector superpo- 
sitions for completeness. 

Conditioned von Neumann measurement. — Consider 
a von Neumann measurement [I], which consists of an 
impulsive interaction Hamiltonian of the form Hj = 
gS(t — to)A <£) p, where A is an observable on the sys- 
tem Hilbert space that we wish to measure and p is 
the transverse momentum on a detector Hilbert space. 
Solving the Schrodinger equation ihdtU = H[U with 
this interaction produces the unitary evolution operator 
U g = exp(g^4 £§) p/ih), which generates translations in 
x by an amount gA due to the canonical commutation 
relations [x,p] = ih. 

Now consider the following experimental procedure. 
First, prepare an arbitrary joint state of the system and 
detector, represented by a density operator psd- Second, 
apply the impulsive interaction U g . Third, measure the 
detector position x or momentum p. Finally, condition 
the detector measurements on an arbitrary generalized 
post-selection on the system, which can always be repre- 
sented by a positive probability operator Pf (23] [24] . 

The conditioned detector averages measured in the lab- 



2 




FIG. 1. (color online) (left) A possible implementation of a 
conditioned polarization measurement similar to [2], where 
the length of a birefringent crystal determines the coupling 
strength g. (right) The weak value Re (as)™ corresponding 
to the Hermite-Gauss detector profiles in Fig. [2] with m = 
(solid, red), m = 1 (dashed, blue), and m = 2 (dot-dashed, 
green), obtained by averaging according to Eq. |TTJ. The 
weak limit g — > is identical for all detectors, as is the 
strong limit g — > oo of a classical conditioned average, but 
the specifics of the transition depend on how the detector de- 
coheres the state. The dotted horizontal line is the eigenvalue 
bound of 1. 
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FIG. 2. (color online) Post-selected detector intensities for 
the initial polarization state \ipi) = (cos(77r/8), sin(77r/8)) 
and final post-selection \ipf) = (1, l)/v2, using the first three 
Hermite-Gauss detector modes with a — 2. Averaging these 
profiles produces weak values according to Eq. ( If I and shown 
in Fig. [T] The dashed line indicates the initial detector inten- 
sity. 



oratory will then have the exact form [24] , 
Tr[(P f ®x)f? SD ] _ 

if) 

Tr[(P f ®p)p SD ] 
Tr[(P f ®l D )pf SD ]' 



Ax) = 
f Tr[(P/ 



M = 



(2a) 
(2b) 



where p' SD — UgPsoU^ is the entangled joint post- 
interaction state at a time i > to- 

As written, Eqs. ([2| show that the joint observables 
Pf <S> x and Pf <8> p are averaged with respect to the fi- 
nal joint state p'g D . However, we can also express these 
averages in terms of the initial joint state by commuting 
the detector observables symmetrically past the evolution 
operators U g to obtain our primary result, 

f (x) =Re(x) w +gRe(A) w , (3a) 
f (p)=Re(p) w . (3b) 

The averages are exactly characterized by the real 
parts of three generalized weak values |22H24j that are 
of the form ([I]), but are on the joint Hilbert space of the 
system and detector, 
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(4a) 
(4b) 
(4c) 



The pre-selection for each weak value is equal to the 
initial joint state psd, while the post-selection is equal 



to the Heisenberg-evolved joint post-selection operator, 
P'sd = Ug(Pf®iD)U g - As noted before, when P f = l s 
there is no post-selection and the weak values Q will re- 
duce to expectation values as a special case. The higher- 
order detector moments are provided in the Supplemen- 
tary Material [25], and all have similar expansions into 
joint weak values. 

Importantly, these relations hold for any coupling 
strength g, any (possibly entangled) initial joint state 
PsDi an d any generalized post-selection Pf; that is, all 
von Neumann detector (conditioned) averages are exactly 
described by generalized weak values. This important re- 
sult seems to have been missed in the existing literature 
due to the fact that the generalized weak values Q can- 
not be written in a form with projective pre- and post- 
selections as defined originally by AAV [T]. Moreover, 
they explicitly include the detector information, so are 
not solely system quantities. 

Reduced state expressions. — If we prepare a product 
initial state psd — Ps ® Pd, where ps (pd) is the ini- 
tial state of the system (detector), then we can exploit 
the product form of the observables to further simplify 
Eqs. pi- Notably, since [A, U g ] = we can express 



Eq. ( 4a I as a weak value only on the system Hilbert space, 

Tr s [PfAp' s 



(Ay 



TrsiPfP's) 



(5) 



where the pre-selection state p' s is the reduced system 
state after the interaction, p' s = Tr£>[p' SD ], and Trg[-] 
(Tr£>[-]) is the partial trace over the system (detector) 
Hilbert space. All detector information has been ab- 
sorbed into an effective preparation of the reduced system 
state p' s . 
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L m (x) 

1 

1-x 
1 - 2x + x 2 /2 
1 - 3a; + 3x 2 /2 - x s /6 



~2 L m (x) 

2 

4-2x 
6 — 6x + x 2 



TABLE I. Laguerre polynomials L m (x) and their derivatives 
for the first few m. These polynomials appear naturally for 
Hermite-Gauss modes in their Wigner distribution (JoJ) , as well 
as the resulting system operations (|10[) and (|12[). 



FIG. 3. (color online) Reduced polarization states corre- 
sponding to the detector responses in Fig. [2] Given an ini- 
tial state ps = (1 + Y2i riOi)/2 with Pauli operators &i and 
measurement of A — 03 with Hermite-Gauss mode m, the 
post-interaction state from Eq. (10 1 is p' s m — (1 + rs&3 + 
L m [{g/o) 2 ] exp(-(p/cr) 2 /2)(ri(Ti + r 2 a 2 ))/2. Bloch sphere 
distortions are shown with the 03 axis aligned vertically; the 
red dot tracks the initial state chosen in Fig. [2] For m > 
any initial state will experience decoherence oscillations and 
pass directly through the 03 axis before partially recohering. 



Since the joint post-interaction state p' SD is necessarily 
entangled by the interaction, the reduced system state 
Ps m @ wm be mixed. However, for sufficiently weak 
coupling one can approximately neglect the interaction in 
([5]) and substitute the initial system state p' s — > ps- The 
detector response ^ will then be linear in g and match 
the original observation of AAV pQ as an approximate 
special case. 

By introducing the Wigner distribution of the de- 
tector state W D (x,p) = ^ fi jdy{x - y/2\p D \x + 
y/2)e wy l h and its Fourier transform Wjj(x, y) — 
$ dpW D {x,p) e ~ ipy/h = (x - y/2\p D \x + y/2), we can 
express the exact reduced system state p' s in a useful 
and compact form, which is our second main result, 



p' s = / dxW D (x,g&d[A})(p s ). 



(6) 



Here ad[A](J3) = AB — BA is the adjoint left-action of 
A as a commutator operation. 



To directly compare the joint weak values Eqs. (4b) 
and ( 4c ) with ([5| , we also express them within the system 
Hilbert space, 



Re (x) w = 



Re(p) 1 



Tr s[P f X{p s )\ 
Tr s [P f V(p s )} 



(7a) 
(7b) 



by introducing the operations X{ps) — T^D[Ug(ps 
(xp D + p D x)/2)Ul\ and T(p s ) = Tr D [U g (p s ® (pp D 



p D p)/2)Ul] that act upon the initial system state. The 
Weyl-ordered operator products that appear allow us to 
use the Fourier transformed Wigner distribution of the 
detector in |6]) to find compact expressions for these op- 
erations, which is our third main result, 



X{ps) = / dxxW D (x,g&d[A])(p s ), 



V{p s ) = ifr 



d z / dxWo(x, z) 



(8a) 

(Ps), (8b) 



z— >g ad [A] 



ihd gzd[A] PS 



Notably, Eq. (8b) allows one to simply obtain the mo- 



mentum response once the functional form of reduced 
system state ^ is known. Generalizations to higher- 
order detector moments are provided in the Supplemen- 
tary Material [25] , 

Hermite-Gauss modes. — To show how our general re- 
sults in Eqs. Q, and ^ can be applied, we 
now consider the Hermite-Gauss modes {|/i m )}, which 
are a widely used complete set of transverse modes nat- 
urally generated in laser cavities that can describe an 
initial zero-mean and collimated detecting beam. The 
Wigner distribution for a Hermite-Gauss mode of order 
rn e {0, 1, 2, ... } has the form [25] , 



W* G (x,p) 
G(x,p) 
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h 2 



(9b) 



where L m is a Laguerre polynomial of order m. The first 
few such polynomials are shown in Table [I] for reference. 

After Fourier-transforming Eq. ^ and integrating ac- 
cording to Eq. ([6]), we obtain a compact expression for 
the exact post-interaction reduced system state for any 
coupling strength and initial detector mode m, which is 
our fourth main result, 



Ps,r. 



-2eC[A] 



s ■ 



(10) 



Notably, a measurement strength parameter e = {g/2a) 2 
naturally appears for all modes along with the Lindblad 



operation C[A] 



-ad [A]/ 2 that decoheres bases or- 



thogonal to the eigenbasis of A [5U 127] ■ Furthermore, 
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the functional form of ( 10 ) is the same as the Wigner 
distribution (|9| up to normalization, but with the func- 
tion G(x,p) replaced by the Lindblad operation — e£L4]. 
Superpositions of modes are considered in the Supple- 
mentary Material [25] . 

Using Eqs. (jToJ, and ^ we obtain the following 
compact results for the exact detector averages for any 
initial Hermite-Gauss detector mode of order m, which 
is our fifth and final main result, 



9 Re (A) w , 



r(p) = 9 



{2a)'- 



■2Im((Ar + A m ). 



(11a) 
(lib) 



Perhaps surprisingly, they are completely parametrized 
by a single generalized system weak value ^ with pre- 
selection equal to the reduced post-interaction system 
state p'g m given in Eq. (10), and one additional weak- 



value-like correction term for the higher mode numbers 
m > 1, 



A m 
M m (ps) = 



Ti s [PfAM m (ps)} 
-2L' m \-2eC[A]\ e^\p s ) 



(12a) 
(12b) 



The first few polynomials —2L' m (x) in Ai m that con- 
tain the derivatives of Laguerre polynomials are shown 
in Table U for reference. 

The appearance of a correction to Im (A) w in Eq. ( lib ) 
further strengthens the observation in [12l [24] that 
Im (A) w pertains solely to the rate of change of the post- 
selection probability and not to the measurement of A 



itself. Indeed, for m = Eqs. (10 1 and (11) correctly re- 



produce the exact Gaussian detector case that we derived 
in more detail using a different method in [24] , 

We stress that these are general results for any system 
observable A. Figs. 

§ and § show the special case 
of an optical application, where A = is a polariza- 
tion observable being measured by a Hermite-Gaussian 
beam. Fig.[2]shows post-interaction spatial intensity pro- 
files for the detector, while Fig. [3] shows the correspond- 
ing reduced polarization states. Fig. [I] shows a possible 
implementation of this example that is analogous to the 
experiment performed in [2], as well as how the general- 
ized weak value ([5]) continuously changes into a classical 
conditioned average as the initial state decoheres. 

Conclusions. — Throughout the controversial history 
of the quantum weak value ([I]), it has been tacitly as- 
sumed that it was a peculiarity specific to the AAV weak 
measurement regime. We have shown in this letter that 
such an assumption has been unwarranted. Indeed, we 
have shown that all (conditioned) averages for any von 
Neumann detector ^ will be completely characterized 
by three generalized weak values Q on the joint Hilbert 
space of the system and detector, which makes such 



weak values a universal feature of von Neumann mea- 
surements. 

We have also shown how to obtain practical and com- 
pact operational expressions for these weak values on 
the system space alone in terms of the reduced post- 
interaction system state Q and two additional opera- 
tions Q. In the process, we have highlighted the prag- 
matic importance of the Fourier transformed Wigner dis- 
tribution of the detector for describing how the detector 
decoheres the system due to the interaction. 

Finally, we have shown that for arbitrary Hermite- 
Gauss modes of a beam-like detector, we obtain simple 
and intuitive operational expressions for the reduced sys- 



tem state ( 10 ) and the (conditioned) detector averages 



(11) that involve the Lindblad decoherence operation. 



The detector averages contain only the real and imag- 
inary parts of a single system weak value ([5]), along with 
a correction ( 12 ) to the imaginary part that appears only 



for the momentum average with higher-order modes. 

Not all observable measurements use such a von Neu- 
mann detector, and not all von Neumann detectors op- 
erate impulsively on the time scales of the system or the 
detector. However, a sufficiently wide class of observable 
measurements use such an impulsive von Neumann pro- 
cedure that the original weak value paper [T] dubbed it 
the "standard measuring procedure," so its universal de- 
scription with generalized weak values is important. We 
also stress that the generalized weak value appears un- 
der reasonable conditions even in the general treatment 
of observable measurements that we developed in detail 
in [23 , a fact which warrants further scrutiny in light of 
the universality shown here. 
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09-0-01417. 
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SUPPLEMENTARY MATERIAL 
Detector Moments 

All detector moments can be determined through con- 
ditioned characteristic functions, 



M Xx ) = 



Tr[(P f ®e iX *)U g pspU}] 
Tr[(P/ ® 1 D ) UgPsoUt] ' 
TT S [P f e^ A X x ( PsD )] 

Tr[(P f <g> 1 D ) U gPsD Ul] ' 
Tr s [P f Vx(psD)} 
TrsiPfP's] ' 



(13a) 



(13b) 



where we have used the Weyl relation 
e iab e —ibp / h e iax ^ anc j nave d e f[ n ed the post- interaction re- 
duced state p' s = r TrD[UgPsDU g '}, as well as the A- 
dependent operations, 

Xx(psd) = \^ D [U g {e^p SD + PsD e^)Ul], (14a) 

Vx(psd) = ^Tr D [U g {e m Psd + Ps D e iXp )Ul], (14b) 

Computing derivatives of the characteristic functions 
produces the conditioned detector moments, 



d" 



f (p n ) 



d(i\y 

gn 

d(i\) n f 



M Xx ) 

(e* Xp ) 



A=0 



A=0 



(15a) 
(15b) 



This procedure is similar in spirit to the full counting 
statistics approach employed in |19j . 

The first two moments are given explicitly by, 



Re(x) w +gRe(A} 1 



(16a) 
(16b) 



f (p) = Re (jp) w , 
f (x 2 ) = Re (x 2 ) w + 2g Re (xA) w + g 2 Re (A 2 ) w , 

(16c) 

f (p 2 )=Re(p 2 ) w , (16d) 

in terms of the Heisenberg evolved joint post-selection 
P' SD = Ug(Pf <8> 1d)U 9 and the joint weak values, 

w = g^MM, (17a) 

^[P'sdPsd] 

{A)W = TtlPkp (A ® i D ) psd] ^ (1?b) 
Tt {P S dPsd] 

<p>* = ^fefa^)^], (17c) 



(x 2 y 

(Ax) 1 
(A 2 ) 1 

(p 2 Y 



TijP'sp (is ® x 2 ) psd] 
^[P'sdPsd] 

Trfe (A ® x) psd] 
^[P'sdPsd] 

Tt[P! sd (A 2 ® 1) psp] 
Tt {P'sd Psd] 

Tr[P SD (i s ®p 2 )psD] 
^[P'sdPsd] 



(17d) 
(17e) 
(17f) 
(17g) 



Detector Wigner Function 

Assuming an initial product state psd 



PS ® Pd, 

we can compute the operations (14) as follows. After 
computing the detector trace in the p-basis and inserting 
two complete x-basis sets, the V\ operation takes the 
form 



P\(psd) 



dpdxdx' 



{x \pd\x) 



2nh 

e -i^{x-x'-gad[A\+h\)fp g \ 

6(x - x' - ,gad[A] -I- hX)(p s ), 
= f dzW D (z,gad[A]-h\)(p s ). 



(18) 



Here we have changed integration variables to z = x — x' 
and y — (x + x')/2, and have noted that Wp(z,y) = 
(z — y/2\po\z + y/2) is the Fourier-transformed Wigner 
function of the detector. 

Performing a similar computation for X\ yields, 



X\(psd) 



dpdxdx' 1 



(e iXx ' + e iXx )(x'\p D \x) 



2irh 2 



(19) 



1 



Tt[P' sdP sd] 



dxdx' ^{e lXx ' +e iXx )(x'\p D \x) 
5{x-x' -gaA[A]){p s ), 
dze iXz W D (z, ga d[A})cos (^ad[A] ) (p s) 



Taking derivatives with respect to (iX) produces the ex- 
pressions in the main text for the first moments. Setting 
A = in either V\{ps) or X\{ps) produces the post- 
interaction reduced system state p' s - 

The operation ad [A] is linear, so any analytic func- 
tion of ad [A] may be defined via its Taylor series in the 
same manner as an analytic function of a matrix. In- 
deed, to more rigorously perform the above derivations 
one can exploit an isomorphism that maps ps into a vec- 
tor and ad [A] into a matrix acting on that vector. After 
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expanding the expressions into the eigenbasis of the ma- 
trix of ad [A] and regularizing any singular functions into 
limits of well-behaved analytic functions, the above com- 
putations can be performed for each eigenvalue, summed 
back into a matrix, and then mapped back into the op- 
erator form shown. For unbounded A then one must 
also carefully track the domains to ensure that the re- 
sulting expressions properly converge, as they should for 
any physically sensible result. 

Hermite-Gauss superpositions 



Notably, the diagonal elements of this sequence are 
the Laguerre polynomials, = ml L m (x 2 ). 

These results can be obtained by using the generat- 
ing function for the Hermite polynomials exp(2xz — 
z 2 ) = X^^Lo H n (x)z n /n\, as well as the identities 
H n (x) — exp(— d2 X )(2x) n and / dxe~ x H n (x)H m (x) — 
6 m , n m\^r2 m . 

Computing the reduced system state p' s using this 
Wigner function yields, 



The Wigner distribution for an arbitrary superposition 
of Hermite-Gauss modes \ip) — ^ m c m \h m ) can be com- 
puted to find (suppressing arguments for compactness), 

W= ^== f fc D™[V2G]e~ G , (20) 



y/ ml n\ 



where G(x,p) = x 2 /2a 2 + 2p 2 a 2 /h 2 and <f>(x,p) — 
tan _1 (— 2pa 2 /hx). The polynomial sequence D™(x) has 
the generating function, 



exp(zz — x(z — z)) = 



j- } D-{x), (21) 



m,n— 



and the explicit form, 



i(m,n) 



k=0 



mini (-l)™- fe 
(m — k)l (n — k)l kl 



(22) 



Ps 



E 



rn.n—O 



Vead[A]\ e^ A \p s ), (23) 



where e = (g/2a) 2 , and C[A] = -ad 2 [i]/2 is the Lind- 
blad decoherence operation. Notably, the functional form 



of the Wigner distribution (20) is still largely preserved 



in Eq. (23 I 



The detector averages can also be computed from this 
Wigner function. The weak value Re (x) w will vanish by 
symmetry; the weak value Re (p) involves the derivative 
*^gad[A]^s> an d, the weak value He(A) w involves the 
state p' s directly. When c m — 1 with the rest of the 
coefficients zero, then these generalizations reduce to the 
results presented in the main text. 



